Generic behaviours in impact fragmentation 

N. SatoiQ, S. Mechkov** and F. Sausset* 

* Laboratoire de Physique Theorique de la Matiere Condensee, Universite 
Pierre et Marie Curie (Paris 6), UMR CNRS 7600 - 4 place Jussieu, 75252 

Paris Cedex 05, France, 

* Laboratoire de Physique Statistique, Ecole Normale Superieure - 24 rue 

Lhomond, 75231 Paris Cedex 05, France 

February 4, 2008 
Abstract 

We present a simple numerical model for investigating the general 
properties of fragmentation. By use of molecular dynamics simulations, 
we study the impact fragmentation of a solid disk of interacting particles 
with a wall. Regardless of the particular form of the interaction potential, 
the fragment size distribution exhibits a power law behaviour with an 
exponent that increases logarithmically with the energy deposited in the 
system, in agreement with experiments. We expect this behaviour to be 
generic in fragmentation phenomena. 

PACS numbers: 46.50.+a, 62.20.Mk, 89.75.Da 
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1 Introduction 

From atomic nuclei to supernovae, including plates and rocks, every cohesive 
system can be broken into fragments, provided that the deposited energy is 
sufficiently large compared to its cohesive energy P~j. The ubiquity and the 
fundamental role of fragmentation are reflected in planetary formation [5] and 
geological phenomena [2j like explosive volcanic eruptions 0], as well as in the 
determination of the phase diagram of nuclear matter EJ. Furthermore, a 
major part of activities in building and pharmaceutical industries makes use of 
particle size reduction processes [HE]. 

In spite of this wide variety of length scales and materials, fragmentation 
processes present some generic features. The most striking of them is the 
frequently observed power law form of the fragment size distribution, whose 
exponent r varies from 1 to about 2.5 according to the data provided by ex- 
periments [H [31 El El |9l [101 HTJ [TT1 (231 [241 [251 EUJ [32l [33l [34l |36] and natural 
phenomena EIHUB]. 

Over the last decades, much effort has been devoted to understanding the 
origin of this scaling behaviour and a number of proposals have been advanced. 
The scattering of the values of r casts some doubts [HI EI] on a self-organized 
critical mechanism proposed by Oddershede et al. [11], but raises the question 
of the existence of universality classes in fragmentation [T21 [T3] ■ Thus, collision 
experiments of atomic nuclei exhibit a percolation behaviour of the fragment 
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size distribution (r ~ 2.2) [6], which is typical of the fragmentation of super- 
critical fluids [2]. Likewise, evidence for criticality has been found in hydrogen 
cluster fragmentation [15] . In contrast, the definition of universality classes in 
fragmentation of macroscopic systems is still an open question. 

As matter stands, fragmentation may depend on several parameters like the 
dimension of the space, the shape, size and material of the fragmenting object, 
as well as on the way it is broken and on the amount of energy deposited in 
it. Simple analytical models based on rate equations [TB], binary fissions [5], or 
sequential fragmentation |17[ 118] although elegant and instructive have limita- 
tions and cannot take into account all the above-mentioned parameters. As an 
example, most of the analytical models predict a single exponent r depending 
only on the dimension of the space, which experiments seem to refute. 

Molecular dynamics calculations have been proven to be a powerful technique 
for addressing the role played by the parameters that may influence fragmen- 
tation processes. Furthermore, simulations allow one to study the dynamics of 
an irreversible and strongly out-of-equilibrium process like fragmentation. Until 
now, mainly two kinds of numerical experiments have been proposed. On the 
one hand, explosive fragmentation of a system is simulated by initially allocating 
a centrifugal velocity to the particles, which interact through a Lennard- Jones 
potential |10 [ 113 ] 119 ] [20], On the other hand, various numerical simulations are 
concerned with the impact fragmentation of brittle solids [T3J [T3] (55] |55J 129] . 
Of particular interest, Kun and Herrmann proposed a sophisticated model that 
takes into account elastic, shear and torque interactions between particles of a 
two-dimensional disk which impacts another disk [12] or a hard plate [25]. They 
suggested that fragmentation of solids occurs as a continuous phase transition 
from a damaged state to a fragmented state when the energy deposited in the 
system increases. 

In the same spirit, we propose a simple model providing a general frame- 
work for investigating the role of the different parameters that may be relevant 
in fragmentation processes. In this paper, we study the impact fragmentation 
of a two-dimensional disc of interacting particles with a wall, as a function of 
the energy deposited in the system. The role played by the materials (interac- 
tion potential) and by the size of the disc (number of particles) are especially 
investigated. 

Note that although fragmenting systems are usually three-dimensional, some 
experiments are performed on platelike objects with a two-dimensional experi- 
mental setup 17, 30 . Furthermore, we are here concerned with general problems 
that may be tackled by two-dimensional simulations. 



2 Model 

To provide a generic frame of reference, the fragmenting system is made up 
of particles which represent mesoscopic grains of materials, without specifying 
their exact nature. These particles interact through a two-body central Lennard- 
Jones type potential: 
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where r^ is the distance between particles i and j, and the two constants, e 
and a, set the energy and length scales respectively. Hence, a is typically the 
diameter of the particles. The two exponent parameters, a and b, control the 
range of attraction and vo is used to set the minimum of the potential to — e. In 
order to mimic the cohesive interaction at a mesoscopic scale, we choose a very 
short-range potential (vq = 107.37, a = 80 and b — 78) with a typical range 
of attraction as short as 0.1 a, in addition to the particle diameter. As we will 
see in the following, our results are not sensitive to the particular choice of the 
parameters vq, a and 6, that is to the range of attraction. 

At first, a disc of a given number of particles is cut in a two-dimensional 
triangular lattice with a lattice spacing corresponding to the minimum of the 
interaction potential given by eq. ([T]). The fragmenting disc is thus a perfect 
crystal (without any quenched disorder) prepared at zero temperature and at 
the corresponding number density of pa 2 ~ 1.16. Next, the disc is rotated by a 
random angle and launched towards a wall, as shown in fig. [T] by assigning to 
each particle a given impact velocity V, perpendicular to the wall. We use the 
magnitude of the impact velocity as a natural control parameter. The particles 
of the disc interact with the wall through the repulsive part of the potential 
given by eq. ([!]), that is vq€(cj / Zi) a , where Zi is the distance between the wall 
and the particle i. We checked that our results do not depend on the precise 
form of this hard-core potential. 

To study the fragmentation process and its outcomes, we perform molecular 
dynamics simulations at constant energy using the Vcrlct algorithm [21 . The 
time step is St = 0.0005 to , where to — \J ea 2 /m is the unit of time and m the 
particle mass, which ensures the conservation of the total energy within 0.01% 
of its average value. In the following, times and velocities are expressed in units 
of tf) and a/t respectively. As we will see below, at t — 500 (10 6 iterations), 
the fragmentation process has already reached a steady state. The main results 
presented in this paper are obtained with a disc of N = 1457 particles, but 
system sizes up to N — 36 289 were used with a view to determine the finite 
size scaling properties of the system. For N — 1457, the potential energy per 
particle is — 2.89e, slightly in excess of the bulk value (— 3e), because of surface 
effects. 

During the fragmentation process, fragments arc identified as being self- 
bound clusters of particles [TU [22] . More precisely, two particles are linked if 
their relative kinetic energy is lower than the absolute value of their interaction 
energy. Note that a definition of clusters based on an arbitrary cutoff distance 
provides the same results when the fragments are far away from each other. 

For a given value of the impact velocity, we calculate the number of frag- 
ments made up of s particles. The fragment size distribution, n(s), is averaged 
over 1000 runs by uniformly sampling the initial random angle of rotation. In 
fragmentation studies, it is usual to estimate the cumulative fragment size dis- 
tribution divided by the fragment size s |llj : 



If n(s) ~ s T , then N(s) exhibits the same, but less noisy, power law behaviour. 
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3 Results 



First we describe qualitatively the collision of the disc with the wall, as the 
impact velocity increases. At very low velocity, a slightly elastic deformation is 
observed at the contact with the wall, followed by the rebound of the disc with- 
out any damage. As the velocity is increased, deformation becomes irreversible 
and fragmentation occurs for V > 0.5. 

Three snapshots of a typical fragmentation event are shown in fig. [T] A zero 
initial angle is chosen in order to clearly visualize the cracks and the compression 
waves materialized by the potential energy of the particles. As can be seen in 
the top snapshot, when the disc strikes the wall, a compression wave starts to 
propagate from the impact zone through the disc. If the impact velocity is large 
enough, the compression wave produces cracks (see the middle snapshot). When 
the cracks coalesce and reach the surface of the disc, fragmentation takes place. 
The part of the initial kinetic energy which is not used to break up the system is 
mainly transformed into the internal and kinetic energies of the fragments. As a 
result of the collision, fragments fly away from each other, as illustrated by the 
bottom snapshot. Note that in spite of the zero initial angle and the symmetry of 
the lattice, irregular cracks give birth to rather rough fragments. To summarize, 
the crack pattern we observe is very much like those found in more sophisticated 
simulations [28, 29 and in two-dimensional experiments [30j H7]. 

It is interesting to note that the speed of sound, c SOU nd, in this material is 
much larger than the impact velocities investigated in this work (0.5 < V < 4.5). 
Indeed, by observing the propagation of the compression waves (see the top 
snapshot in fig. [I]), we roughly estimate c SOU nd — 100. Furthermore, a calculation 
of the speed of sound in the harmonic approximation for a one-dimensional 
chain of particles interacting through the potential v(r) gives the same order of 
magnitude: c sounc i = 78. 

In order to better understand the nature of this fragmentation process, we 
calculate the mean size of the first and second largest fragments, denoted S max i 
and <S , maX 2 respectively. As with percolation theory, we also compute the average 
of the second moment of the fragment size distribution, which is related to the 
mean fragment size |31j : 

rn 2 = ^ s 2 n(s), (3) 

s 

where the sum excludes the largest fragment. 

These three quantities are plotted in fig. [2] as a function of the impact ve- 
locity, at t = 100 and at the end of the simulation at t = 500. We clearly see 
that S'maxi, Smax2 and ?7i2 evolve little between these two times, showing that 
the fragmentation process has reached a steady state. Indeed, the fragment 
size distribution is rapidly stabilized after the impact. Although the fragments 
are heated during the process, their internal energies are low and allow the 
evaporation of only a few particles. 

As expected, S'maxi decreases with V, reflecting the increasing violence of the 
impact. On the other hand, both S , ma x2 and mi have a maximum at the same 
threshold velocity, Vt — 1.15, as in the simulations performed by Kun, Herrmann 
and their collaborators [TH H5J . These authors interpreted this behaviour as a 
sign of criticality in fragmentation. We note that the value of Vt decreases 
slightly with TV (for instance, Vt — 1 for N = 36 289). As is shown in fig. [2] the 



4 




Figure 1: (colour online) Snapshots of the fragmentation of a disc of N — 36289 
particles launched with a zero initial angle and an impact velocity V = 2, at 
t = 1 (top), t = 10 (middle) and t — 100 (bottom). Particles are coloured 
according to their potential energy from — 6e (blue) to (red). 
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Figure 2: The first and second largest fragments and the second moment as a 
function of V for a disc of N = 1457 particles at time t — 100 (open symbols) 
and at t = 500 (solid symbols) . The full lines correspond to power law fits (see 
text). The threshold velocity is indicated by an arrow. 



decrease of S ms ,xi, <Smax2 and is well fitted by a power law ~ V~ J , with an 
exponent 7 equal to 1.2, 1.6 and 2.6 respectively. As can be remarked, in contrast 
with a critical behaviour, the threshold velocity does not appear in this power 
law form. Moreover, the fit is very good not only in the vicinity of Vt, but also 
at velocities up to 3.5. It is worth noting that the power law behaviour of S max i 
has been reported in experimental works dealing with cubic ice (7 = 1.8) [32] 
and aluminum spheres (7 = 6.6) [33 . Moreover, by fragmenting various rock 
types, Matsui and his coworkers [23] showed that the exponent seems to depend 
mainly on the shape of the system: with 7 ~ 1.4 for spherical samples and 7 ~ 3 
for cubic ones. 

For V > V t , the fragment size distribution N(s) follows a power law with 
an exponent t (see fig. |3|. Of course, the properties and the number of small 
fragments (s < 10) may depend strongly on the particular features of the sys- 
tem [T2"| 158], However, we are mainly concerned with the large size range for 
which generic behaviours are expected. 

To test the dependence of our results on the potential interaction, we carry 
out simulations with two other potentials having the same form given by eq. (JT|) , 
but with a longer attraction range of a (vq = 16.29, a = 13 and b = 11) and 2a 
(vq = 4, a = 12 and 6 = 6), the latter being a standard Lennard- Jones potential 
(see the inset of fig. [3| . The fragment size distributions corresponding to these 
three potentials are plotted in fig. [3] Strikingly, they share the same slope, and 
so the same value of t, suggesting the robustness of our results. Differences 
occur in the small size region and for the largest fragment, the size of which 
increases with the range of the potential, as well as the cohesion of the disc. 
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Figure 3: Fragment size distributions at V = 2 for a disc of N = 1457 particles 
and for three different interaction potentials v(r) plotted in the inset. The 
attraction range is given in the legend. 



Furthermore, this result implies that fragmentation may be, to a certain extent, 
independent of the particular material as suggested by Oddershede et al. . 

As mentioned above, the fragment size distribution exhibits a power law 
behaviour not only at Vt , but for all the values of the control parameter greater 
than V t . What is more, we find the exponent r to grow logarithmically with the 
impact velocity (see fig. [4| as 

T = a\nV + f3 1 (4) 

with a = 0.70 ± 0.01 and (3 — 1.00 ± 0.01. Because of the decrease of 5 max i 
with increasing V, the power law region narrows, which enlarges the errors bars. 
Nevertheless, the logarithmic fit is quite sensible for Vt < V < 3. 

The growth of the exponent r with the energy imparted to the fragmenting 
system is clearly observed in various materials of different shapes, like glass and 
plaster plates [3D] , ceramic tubes [33] , gypsum disks [IB] and glass rods [TD] [55] . 
However, as far as we know, only a few experimental papers are concerned with 
the exact dependence of r on the energy. Thus, by performing impact experi- 
ments on rocks, Matsui and his collaborators [24] found the logarithmic relation 
given by eq. Q with a ~ 0.65, regardless of the shape and rock type investigated 
(granite, basalt, tuff). Furthermore, Moukarzel et al. have very recently studied 
the fragmentation of liquid (glycerol and water) droplets by a pressurized-gas 
blow [23] . As can be seen in fig. 6 of ref. [23] , r has approximately a logarithmic 
behaviour as a function of the jet pressure, which is related to the imparted en- 
ergy. The increase of r with energy reported in other simulations of impact [55] 
and explosive [TU] fragmentation, seems to be in agreement with eq. Q , at least 
at low impact energy. 



7 




Figure 4: Exponent r of the fragment size distribution as a function of the 
impact velocity for a disc of N — 1457 particles. The full line corresponds to a 
logarithmic fit (see eq. Q). 



In contrast, by performing the same kinds of numerical experiments, other 
authors reported a virtually constant value of r with increasing impact en- 
ergy p2J [TH [201 EE], corresponding to a zero value of a. Also, as mentioned 
above, simple analytical [HI [T71 [TH] and numerical models [27J [55] predict a value 
of r which depends only on the dimension of the space. 

Considering these experimental and theoretical results, we suggest that the 
logarithmic dependence of the exponent r on the impact energy (or velocity) 
given by eq. Q is generic of fragmentation processes. A question which natu- 
rally arises is on what parameters the coefficient a depends. 

We now turn to the finite-size scaling properties of the fragment size distri- 
bution. It follows from the conservation of the total number of particles that: 

N 

i\T = 5>n(s). (5) 

s=l 

Suppose that n(s) ~ f(N) s~ T , where / is a function of N to be determined. 
By taking the continuous limit (large N), we have 

N ~ /(AO J s 1 ^ ds ~ f(N)N 2 -\ (6) 

as long as r < 2. Hence, f(N) ~ jV r-1 , and the fragment size distribution 
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Figure 5: (colour online) Scaling collapse of the fragment size distribution at 
V = 2. From right to left N = 89, 515, 1 457, 5 813, 10 181, 36289. In the inset: 
Smaxi and 77i2 are fitted by N (solid lines) and N 2 (broken lines) respectively. 



scales as: 

n(s) ~ N(s) ~ 1(A)" 7 '. ( 7) 

The fragment size distribution N(s) defined by eq. ^ is plotted in fig. [s] for 
various sizes N of the disc as a function of s/N. Considering the important 
surface effects at small N, the data collapse is good, except for the small frag- 
ment sizes, which do not follow a power law with exponent r. We find that the 
quality of the data collapse is better at higher impact velocities. 

The validity of eq. Q is confirmed by the scaling behaviour of the second 
moment of the fragment size distribution given by eq. p}. Using eq. Q, we 
obtain (with r < 3): 

/•N 

m 2 ~N T - 1 J s 2 - T ds~N 2 . (8) 

As illustrated in the inset of fig. [5] the scaling behaviour predicted by eq. ^ 
is in very good agreement with the simulation data. Furthermore, the size of 
the largest fragment scales nicely with the total number of particles, that is 
Smaxi ~ N. It must be emphasized that the scaling behaviour of S'maxi and mi 
is independent of t, and is then observed for all the velocities larger than V t 
investigated in this work. 
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4 Discussion 



To summarize, this simple model provides a general framework for investigating 
the various parameters that may play a role in fragmentation processes. As we 
have seen, our results do not depend on the range of attraction of the interaction 
potential, which suggests a certain universality in fragmentation of solid bodies. 
Qualitatively, the propagation of cracks and the crack pattern obtained in our 
simulations are akin to the ones observed experimentally in platelike objects. 

The power law behaviour of the fragment size distribution and the maximum 
of its second moment seem to provide evidence that fragmentation is a contin- 
uous phase transition, similar to a percolation transition, as suggested by Kun 
and Herrmann [TJ]. However, in contrast to continuous phase transitions, the 
power law is observed, not only at the threshold velocity, but also at all veloci- 
ties above Vj- Moreover, the scaling behaviour of the fragment size distribution 
is simply inferred from the conservation of the number of particles and does not 
reveal a critical behaviour. We note that these properties are very similar to 
the avalanche phenomena associated with hysteresis loops, in which Sethna and 
his collaborators |37j found scaling and power law behaviours in a large region 
near the critical point. However, this analogy has to be confirmed by thorough 
analysis. 

Furthermore, we find that the exponent of the fragment size distribution 
increases as a logarithmic function of the imparted energy, in agreement with 
experiments on rocks and on liquid droplets. This result casts some doubts 
upon the existence of universality classes in fragmentation, in the strict sense of 
continuous phase transitions, but it sheds light on the scattering of the values 
of t measured experimentally. Besides, we expect this logarithmic relation to 
be generic in fragmentation phenomena. 



We thank H.J. Herrmann for helpful conversations, P. Viot and H. Hietala 
for critically reading the manuscript. 
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